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Abstract 

In this paper, a non-Lyapunov novel approach is proposed to estimate the 
unknown parameters and orders together for non-commensurate and hyper 
fractional chaotic systems based on cuckoo search oriented statistically the 
differential evolution (CSODE). Firstly, a novel Gao's mathematical model is 
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put and analysed in three sub-models, not only for the unknown orders and 
parameters' identification but also for systems' reconstruction of fractional 
chaos systems with time-delays or not. Then the problems of fractional- 
order chaos' identification are converted into a multiple modal non-negative 
functions' minimization through a proper translation, which takes fractional- 
orders and parameters as its particular independent variables. And the ob- 
jective is to find best combinations of fractional-orders and systematic pa- 
rameters of fractional order chaotic systems as special independent variables 
such that the objective function is minimized. Simulations are done to es- 
timate a series of non-commensurate and hyper fractional chaotic systems 
with the new approaches based on CSODE, the cuckoo search and differen- 
tial evolution respectively. The experiments' results show that the proposed 
identification mechanism based on CSODE for fractional-orders and param- 
eters is a successful methods for fractional-order chaotic systems, with the 
advantages of high precision and robustness. 

Keywords: Non-commensurate and hyper fractional-order chaotic system, 

unknown fractional orders, unknown systematic parameters, cuckoo search 

oriented statistically the differential evolution, non-negative special 

functions' minimization 

PACS: 05.45.Pq, 02.60.Pn 
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1. Introduction 

The applications of fractional differential equations began to appeal to 
related scientists recently[l-27] in following areas, bifurcation, hyperchaos, 
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proper and improper fractional-order chaos systems and chaos synchronization[l- 
33]. 

However, there are some systematic parameters and orders are unknown 
for the fractional-order chaos systems in controlling and synchronization. It 
is difficult to identify the parameters in the fractional-order chaotic systems 
with unknown parameters. Hitherto, there have been at two main approaches 
in parameters' identification for fractional-order chaos systems. 

• Lyapunov way. There have been few results on parameter estimation 
method of fractional-order chaotic systems based on chaos synchronization [34] 
and methods for parameter estimation of uncertain fractional order 
complex networks[35]. However, the design of controller and the updat- 
ing law of parameter identification is still a tough task with technique 

and sensitively depends on the considered systems. 

• non-Lyapunov way via artificial intelligence methods For examples, 
such as differential evolution[7] and particle swarm optimization [9]. In 
which the commensurate fractional order chaos systems and simplest 
case with one unknown order for normal fractional-order chaos sys- 
tems are discussed. However, to the best of our knowledge, little work 
in non-Lyapunov way has been done to the parameters and orders esti- 
mation of non-commensurate and hyper fractional-order chaos systems. 
And there are no general mathematical model has been purposed for 
all these kinds of identification. 

We consider the following fractional-order chaos system with time delays. 

a D«Y(t) = f(Y(t),Y (t),9,T) (1) 
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where Y(t) = (yi(t),y 2 (t), ...,y n (t)) T G 3? n denotes the state vector. 9 = 
(9\, 6> 2 , 9 n ) T denotes the original parameters, r — (tl, t 2 , r n ) is the time 
delay, g = (qi, q 2 , q n ), (0 < ^ < 1, i — 1, 2, n) is the fractional derivative 
orders. 

f\(Y(t),Y (t),e,r) = (/l'/2>->/n)l(y(t),y (t),ff,T) 

Normally the function / is known. And the 0, g, r are unknown , then 
the = 6*2, 9 n , qi, q 2 , q n , Tl, 7"2, Tn) will be the parameters to be 
estimated. 

Then a correspondent system are constructed as following. 

a D q t Y(t) = f(Y(t),Y (t),9,f) (2) 

where Y(t),9,q,f are the correspondent variables to those in equation (1), 
and function / are the same. The two systems (1) (2) have the same initial 
condition Y (t). 

Then the objective is obtained as following, 

6* = arg min F = arg min 1 1 Y (t) — Y (t) (3) 

t 2 

When some the fractional chaotic differential equations / = (/i, / 2 , /„) 
are unknown, how to identify the fractional system? That is, 

(/i, f2, fn)* — arg min F (4) 

Ul,J2, — Jn) 

Now the problem of parameters estimation (3) become another much more 
complicated question, fractional-order chaos reconstruction problem [36], to 
find the forms of fractional order equations as in (4). In Reference [36] a novel 
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non-Lyapunov reconstruction method based on a novel united mathematical 
model was proposed to reconstruct the unknown equations (/i, fi, f n )- 

When it comes to the system (1) a D q t Y(t) = f{Y(t), Y (t), 9) with neither 
q nor 8 are known, the united model are not effective. For the united math- 
ematical model[36], to be identified is only (f 1: f 2 , /„) instead of q. That 
is a DtY(t) of the equation (2) are not included. Actually, if the q are taken 
into consideration in the united model (4), then the basic parameters' setting 
to be reconstructed in Reference [36] will be basic set {x, +, — } with ex- 
tra } etc., and the input variables {x±,X2, ■ ■■ 1 x n } with {Y} extra and 
etc. Although for the candidates "programs" in Reference [36] the maximum 
depth of tree is 6, considering the maximum number of nodes per tree is in- 
finite, there will be infinite illegal candidates will be generated. Then, in one 
hand, the most time-consuming thing for the novel united model (4) is to kill 
these illegal individuals from the legal individuals. However, these defaults 
are not solved in Reference [36]. In the other hand, as q G D\ in unknown, 
it is really difficult to generate an individual with {x, -=-,+,—,=, Df}, in 
neither illegal nor legal cases. And up till now, there is no existing way to 
resolve these defaults. And we can conclude from simulations [36] that the 
proposed method are much more efficient for the systems with coefficients in 
(fi, f2, fn) as integers than as improper fractions. 

Therefore, to estimate the q of the equation (2) with unknown system- 
atic parameters 9 is still a question to be solved for parameters and orders 
estimation of non-commensurate and hyper fractional-order chaos systems. 

And Cuckoo search (CS) is an relatively new and robust optimization 
algorithm[37, 38], inspired by the obligate brood parasitism of some cuckoo 
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species by laying their eggs in the nests of other host birds (of other species). 
The searching performance is mainly based on the Levy flights mathematically[37- 
39], which essentially provide a random walk while their random steps are 
drawn from a Levy distribution for large steps[37-39]. However, in of CS 
evolution, the Levy flights in each main iteration are used twice. It has two 
results, the CS' searching performance become a little strong but the re- 
dundant evaluations for the objective function are generated too. Therefore, 
some more efforts are need to improve the performance of CS. 

To the best of authors' knowledge, there are no methods in non-Lyapunov 
way for non-commensurate and hyper fractional order chaotic systems' pa- 
rameters and orders estimation so far. The objective of this work is to present 
a novel simple but effective approach to estimate the non-commensurate and 
hyper fractional order chaotic systems in a non-Lyapunov way. And the il- 
lustrative reconstruction simulations in different chaos systems system are 
discussed respectively. 

The rest is organized as follows. In Section 2, a general mathematical 
model not only for fractional chaos parameters identification but also for re- 
construction in non-Lyapunov way are newly proposed and analyzed in three 
sub-models A, B and C. And a simple review was given on non-Lyapunove pa- 
rameters estimation methods for fractional-order and normal chaos systems. 
In section 3, a novel methods with proposed united model based on Cuckoo 
search oriented by differential evolution statistically (CSODE) is proposed. 
And simulations are done to a series of different non-commensurate and hy- 
per fractional order chaotic systems by a novel method based on CSODE in 
Section 4. Conclusions are summarized briefly in Section 5. 



6 



2. Gao's mathematical model for fractional chaos reconstruction 
and orders estimation in non-Lyapunov way 

In this section, a general mathematical model for fractional chaos param- 
eters identification in non-Lyapunov way is proposed. A detail explanation 
for the general mathematical model will given as following subsections in 
three aspects, sub-model A, B, and C. 

2.1. Gao's mathematical model 

Now we consider the general forms of fractional order chaos systems (1). 
To make the system (1) more clear, we take its equivalent form as following 
system (5). 



aDfy\{t) = fi(t, yi(t), yi{t - T),y 2 (t),y 2 (t - r), ...,y n (t),y n (t - r)) 
a Dfy 2 (t) = f 2 (t, y^t), y^t - r), y 2 (t), y 2 (t - r), y n (t),y n (t - r)) 

a D q t n y n (t) = f n (t, yi(t), y x (t - r) , y 2 (t) , y 2 (t - r), y n (t),y n (t - r)) 
L = (y 1 ,y 2 ,...,y n ) 

(5) 

And a correspondent system (6) is constructed as following. 

a Df 1 y 1 (t) = f 1 (t,y 1 (t),y 1 (t-f),Ut),Ht-r),...,Ut),Ut-r)) 
a Df 2 y 2 (t) = f 2 (t, yi(t), yi(t - f), y 2 (t), y 2 (t - f), y n (t),y n (t - f )) 

*D q t n y n {t) = j n (t, y 1 (t),y 1 (t - f), y 2 (t), y 2 (t - f), ...,y n (t),y n (t - ?)) 
L = (y u y 2 ,...,y n ) 

(6) 
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To have simple forms, we take a = 0. 

Then the novel objective function (fitness) equation (7) in this paper come 
into being from equations (18) and (19) as below. 

T-h 



F= £ \\ L 



t=0-h 



(7) 



Now a novel Gao's mathematical model for fractional chaos reconstruction 
comes into being as Figure 1 shows, where functions / = (/i, / 2 , /„), 
/ = (/1,/2,-,/n), fractional orders q = (q 1: q 2 , q n ), q = (gi, q 2 , •-, q n ), 
time-delays r = (ri,r 2 , ...,r n ), r = (fi,T2, •••,T n ), systematic parameters = 
(6>i, 6*2, 0„), 9 = (0~i, 2 , -, 0n) respectively. 
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Figure 1: Gao's mathematical model for fractional chaos reconstruction 



A detail explanation for the general mathematical model will given as 
following subsections in three aspects, sub-model A, B, and C. 

And the objective function (7) to be optimized can also be any kind of 
equations (8), (9), (10), (11) to be minimized by artificial intelligence methods 
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as following. 



A? 



1 N ii 

g= nJ2\\ l *- 1 



=1 



N 



H = ^ 1 1 Li — Li 
1 N ii 



(8) 



(9) 



(10) 



;n) 



2.2. Mathematical sub-model A. 

It should be noticed here that the independent variables in function 
(7) in the general model in Figure 1 are not always the parameters and 
fractional orders. They can be the special variables, for instance, as func- 
tions / = (fi, f 2 , f n ), fractional orders q = (qi, q 2 , q n ) and time delays 
t = (fi,f 2 , ...,r n ). 

And for the sub-model A, that is 



arg min F 



JA,f 

It can also be written as following. 

((A, h, fn), (qi, 92, -, Qn), (fi, f 2 , f„)) = 



(12) 



arg min F 

((/l,/2,-,/n),(91>92,.-->9n),(n,f 2 ,...,f„)) 



There exist several definitions of fractional derivatives. Among these, the 
Griinwald-Letnikov (G-L), Riemann-Liouville (R-L) and the Caputo frac- 
tional derivatives are the commonly used[40-45]. And G-L, R-L and Caputo 
fractional derivatives are equivalent under some conditions [46]. 
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The continuous integro-differential operator [47, 48] is used, and we con- 
sider the continuous function f(t). The G-L fractional derivatives are defined 
as following. 



h 



a D<f(t) = Jim 1 g (-1)' |j J / (f - (13) 

where [x] means the integer part of x, a,t are the bounds of operation for 
a D q J(t), q e 3*. 

We take ideas of a numerical solution method [47, 48] obtained by the 
relationship (13) derived from the G-L definition to resolve system. That is, 

{k _ Y) D*J(t) ~ I £ (-1)^ * j / (f fc _,-) 

where L m is the memory length, t fc = fc/i, h is the time step of calculation and 
( — iy I ) are binomial coefficients cf \ (j = 0, 1, ...,). When for numerical 

\ J / 

computation, the following are used, 

Then in general, for the simplest case (14) of equation (5) as following. 

DU(y(t)) = f(t,y(t),y(t-r)) (14) 

Let ykr — y(tk — T(tk)) -It can have the approximate value as equation 
(15), when it used for calculating. 

r/h-[r/h] l + [r/h]-r/h 

ykr= Vk+l+lr/h] ^ 1" yk+[r/h] ^ (15) 
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And let y k = y(t k ) , then equation (14) can be expressed as 



k 

,(«). 



j=0 



y (4) = / (4, ?/ (4) , y (t k - r(t k ))) h q -J2 cfy (t k -j) 

j=v 



yk = f {t k , yk, Vkr) h q - cfvk-j (16) 

j=v 

where v in above is defined as 



v = 



K ~h~i IT 

or v — 1 for all k. 

Equation (16) is a implicit nonlinear equation respect to y k . Now we can 
construct an iteration algorithm to solve y k as following (17). 



J=V 



where / is the iteration number. When y k — y k < 5 ( normally the 
given error 5 < 1CT 6 ), we consider y k = y k +1 ^ be the solution of the simplest 
equation (14). And if the derivative of the / exist and |/| < M (M is 
constant), and h q M < 1, then the iteration (17) converges to a constant as 
long as the calculus step h is smaller enough. 

With the ideas from iteration of equation (17), the systems (5) and (6) 
are solved as following (18) and (19) respectively, 
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X 



(1+1) 



i k 



fi (tk-1, tk, 



T (0 T (0 T (0 T (0 T (0 T (0 T (0 T (0 ) /,« 

x l fc ' Xl fc,r' k ' fc,r' x * fe-1' x * fc-l,r' x ™ fc-1' x ™ fc-l,T / 



fe 

E 



(18) 



x « fe — 

f ft* i t* 5i (0 5! (,) 5- i (,) 5- i (0 5- (0 5- (0 5 (0 5 (,) ^ h* 

Jl l t-fc-l, t-fc, J'l fc 5-^1 ' fe ' fc,f ' Jjl fe-1' X * fe-l,f •••' Xn fe-1' X ™ fc-l,f J " 

fe ^ / ~\ 

Cj %i k—j 

j=v 

(19) 

And if r = (0, 0, ...,0), then the systems (5) can be solved as equation 
(20) [36]. 

fe 

#i i ^ — /i f^l ^i-l fe\ fc-D •••) x n U-l) ^ ~ C< f k-j (20) 

j=v 

To our best of knowledge, there is no work have been done to reconstruct 
the fractional chaos systems under condition that both /, q and r are un- 
known in sub-model A as equation (12) neither for time-delays free nor with 
time-delays chaos systems . 

2.3. Mathematical sub-model B. 

In this sub-model, / are unknown but r and q are definite. Then to be 
estimated is only the fractional differential equations / , that is 

(/) =argminF (21) 
V ) (/) 
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It is should be noticed that there is few method for reconstruction for 
fractional order chaos systems [36] so far. 

However, there are a few results for normal chaos systems, as the spe- 
cial cases of fractional chaos systems. For reconstruction of / = (/i, f 2 , f n ) 
with the non-Lyapunove methods, they are mainly from symbolic regression 
through genetic programming (GP) [49-51], and some evolutionary algorithms [52- 
62]. 

Considering mathematical sub-model A, we have to say it is really difficult 
to use the ideas in mathematical sub- model B. Let the input variables are 
taken as x, y, z and the basic operators set used be {+, — , x , D\ , =}, where 
fractional order p e [0, 1] is uncertain. Now we consider easiest cases that the 
fractional order differential equation Dfx = fi(x, y, z) is unknown. Then we 
will see the individuals as following with the ideas similar to methods for the 
normal chaos of sub- model B. Figure 2(a) and Figure 2(b) are the normal and 
correct candidate individuals only for the right part f\ of dx/dt = fi(x,y,z) 
of the normal chaos systems. 

However, when it comes to fractional order chaos system, the whole frac- 
tional order differential equations should be taken into accounts, that is 
Dfx = fi(x,y,z) with q = p e [0,1] uncertain and f 1 unknown. Figure 
2(c) shows a correct candidate. And when the evolutions (crossover, mu- 
tation and selection) go on, there are some wrong and illegal candidates 
generated as Figures 2(d) , 2(e), 2(f) show. Figure 2(d) is a wrong candidate 
with Df(x — y). Figure 2(e) is a wrong candidate with Df(x — y) and D^ 2 (y). 
Here it should be noticed that random p±,p2 € [0, 1]. Figure 2(f) is a wrong 
candidate with not only D^(x — y), Df 2 (z) and Df\z) but also extra {=}. 
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Here it should be noticed that random pi,P2,Pz G [0, 1]. 



a a 

sin In 2/ x 



A 



Z X z 



sin x 

4- A A 

x y z x y 

(a) sin x + yz (b) sin(x — y) ■ laz — (y + xz) 

-A. 

A. 7 A 

Df x - y x 

4- / V A / \ 

x y z x y x z 

(c)Dfx = yz (d) D p t {x-y) = y + xz 



X 



£>f + - y - 

/v A a 

y x x 2/ 4- Df 

/V A A 4- 

x y / Df 2 L>f y I 

4- 4 

x z 

(e) Df 1 (x — y) = y + zD P2 x (f) Illegal candidatures 

Figure 2: Some examples of the tree structures in GP evolutions 

So long as the evolutions (crossover, mutation and selection) go on, these 
wrong candidate are inevitably existing, although in the genetic program- 
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ming, the tree depth is set to be limited with unlimited leaves. And these 
kind of wrong individuals will become heavy burden for both the genetic 
evolution and resolving the fractional order differential equations. 

Thus it is not suitable to use the methods based on GP only to reconstruct 
the fractional order chaos system neither fractional order q nor equations fi 
are unknown. However, if it is only considering the unknown equations fi 
with definite certain fractional order q , these methods will be impressive and 
efficient as in Reference [36]. 

2.4- Mathematical sub-model C. 

In this sub- model C, (gi 7^ q2 7^ ... 7^ q n ) , systematic parameters and 
time delays (9 1 , 9 2 , 9 n ), f = (fi, f 2 , f n ) G (fi, f 2 , f n ) are unknown for 
non-commensurate and hyper fractional order chaos system. 

There some estimation methods have been purposed to identify the un- 
known parameters and orders for commensurate fractional order chaotic sys- 
tems. However, to our best of knowledge, no such reconstruction methods 
have been done for non-commensurate and hyper fractional order chaos sys- 
tem, it is necessary to resolve the following (22) in non-Lyapuno+v way. 



And there exist basic hypotheses in traditional non-Lyapunov estima- 
tion methods for fractional order systems[7, 9, 28]. That is, the param- 
eters and fractional orders are partially known or the known data series 
coincide with definite forms / = (/1, / 2 , f n ) of fractional chaotic dif- 
ferential equations except some uncertain parameters and fractional orders 
= (#!, #2, 9 n , qi, q 2 , q n ). 




(22) 



{q,e,f) 
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This is the basic difference between sub-model A,B and C. And for the 
case when some the chaotic differential equations / = (/i, f 2 , ■ f n ) are 
unknown, there are some chaos system reconstruction methods. Then the 
cases [63-75] can be thought as special cases of chaos reconstruction, when the 
exact forms of chaotic differential equations / = (/i, f 2 , f n ) are available 
but some parameters unknown. 

2.4-1 ■ Parameters estimation for fractional order chaos systems 

We take the fractional order Lorenz system (23) [3, 8, 24] for example, 
which is generalized from the first canonical chaotic attractor found in 1963, 
Lorenz system [76]. 



uPfx = a-(y-x); 

— m 

(23) 



a D q t 2 y = 7 • x - x ■ z - y; 



a Dfz = x-y -b- z. 
% L = (x,y,z) 

where qi, q 2 , q% are the fractional orders. When (g 1; q 2 , (fe) = (0.993, 0.993, 0.993), 
a = 10,7 = 28,6 = 8/3, a = 0, intimal point (0.1,0.1,0.1) system (23) is 
chaotic. Generally when the dimension 

= 9i + ?2 + <?3 > 2.91 

for fractional system (23) is chaotic[3, 8, 24]. 

The form of function (12) can also be as following: 

T-h 2 

p 2 =F(a, 1 ,b,q 1 ,q 2 ,q 3 )= ^ ||Z - (24) 

t=0-h 

It is noticed that the objective function (24) can be any forms of corre- 
spond equations (8), (9), (10), (11). 
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Then the problems of estimation of parameters for chaotic system are 
transformed into that of nonlinear function optimization (24). And the 
smaller p 2 is, the better combinations of parameter (a, 7, b, q\, q 2 , (73) is. The 
independent variables of these functions are 9 = (a, 7, b, q±, q 2 , 93). 

And considering the fractional system is very complicated, to simplify the 
problems, it is reported unknown q = qi = q 2 = 53,(7,7,6 or case of a, 7, b 
are known and only one qi are unknown for the similar fractional order chaos 
systems, such as fractional order Lii system[16, 77] fractional order Chen 
system[27, 78] fractional Lorenz system[3, 8, 24], discussed in Ref. [7, 9]. 
The above is the basic idea for the recently proposed methods for fractional 
chaos system[7, 9]. 

However, the case q\ 7^ q 2 7^ (fe are not included in the above non- 
Lyapunov ideas or not fully discussed either for non-commensurate fractional 
chaos systems. 

2.5. The main differences between sub-models A, B and C 

Equation (12) is the crucial turning point that changing from the pa- 
rameters estimation into functions reconstruction and orders estimation, in 
other words, both fractional order estimation and fractional chaos systems' 
reconstruction. 

It can be concluded that the parameters' estimation of fractional order 
chaos system[7, 9] is a special case of fractional order chaos reconstruction 
here as (12). In their researches, the forms of the fractional order differen- 
tial equations f 2 , f n ) are known but some parameters 9 2 , 6 n ) of 
these equations are unknown, and only one the fractional order and some of 
these systematic parameters (0i,0 2 , ■■■,0 n ) are estimated[7, 9]. 
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And further, the parameters estimation cases that all (fi, / 2 , f n ) are 
known but parameters (#1, #2, n ) of these equations are unknown, and the 
reconstruction case that some of {fi, f'2, fn) are unknown, in Section 2.3 
for the normal chaos system, are the special cases of fractional order chaos 
systems' reconstruction (12). 

However, it should be emphasized here that , for reconstruction the novel 
general mathematical model (12) for fractional chaos parameters identifica- 
tion in non-Lyapunov way, with uncertain different fractional order q, that 
is Ql ± q 2 ^ ... ^ q n e {x, +, -, Df,Dl 2 , Df n }, it is really difficulty to 
generate a proper candidate from this basic set as shown in Figure 2. Then, 
it is not easy to reconstruct the fractional order differential equations and 
identify the fractional orders together. And only the simplest case that with 
definite q = qi = q 2 = ... = q n are discussed [36]. 

3. Cuckoo search oriented statistically by differential evolution 

3.1. Cuckoo Search 

Cuckoo search (CS) is an optimization algorithm[37, 38], inspired by the 
obligate brood parasitism of some cuckoo species by laying their eggs in the 
nests of other host birds (of other species). And some host birds can engage 
direct conflict with the intruding cuckoos. 

CS is based on three idealized rules: 

• Only one egg is laid, and is dumped into a randomly chosen nest by 
each cuckoo at time t; 

• The best nests with high quality of eggs (candidate solutions) will be 
copied to the next generation directly; 
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• The number of available host nests is fixed, and an alien egg will be 
discovered by a host bird with probability p a G [0,1]. If so, the host 
can either throw the egg away or abandon the nest so as to build a 
completely new nest in a new location. 

A Levy flight is performed for cuckoo i when a new candidates x^ t+1 ^ is 
generated [3 7-39], 

x f +1 > = x f + a © Levy(A) (25) 

where a > is the step size which should be related to the scales of the 
problem of interest. Normally, a = 0(1). The product © means entry-wise 
multiplications. Levy flights essentially provide a random walk while their 
random steps are drawn from a Levy distribution for large steps 

Levy ~ u = t" A , (1 < A < 3) 

which has an infinite variance with an infinite mean, essentially form a ran- 
dom walk process obeying a power-law step-length distribution with a heavy 
tail [37-39]. 

Based on above rules and ideas, the basic steps of the CS can be sum- 
marised as following pseudo code Algorithm 1. 

It should be noticed that in each iteration of Algorithm 1, there are two 
rounds of evaluation the fitness, ones is after getting a cuckoo by Levy flights, 
the other is after abandon the worse nests with probability p a and building 
the new nest at the new locations. It is also showed in the original Matlab 
code in Reference [38]. 

This might be the reason that CS is efficient. Because CS use Levy 
flights twice and evaluate the candidates twice in one generation. However, 
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Algorithm 1 The basic steps of the cuckoo search 
1: Basic parameters' setting Objective function f(x),x = (x±, Xd), 

initial population of n host nests Xi(i = 1,2, ...,n), boundaries for each 

dimension Xi = (x^i, ...,Xi t d) and etc. 
2: while Termination condition is not satisfied do 
3: Get a cuckoo (i, i = 1,2, n) randomly by Levy flights (25). 
4: Evaluate its fitness F^, 
5: If (F t > Fj) 

6: Replace j by the new solution 
7: End 

8: Abandon a fraction (p a ) of worse nests. 

9: Build new ones at new locations via Levy flights (25) (i, i = 1, 2, n) 

10: Keep the best solutions (or nests with quality solutions). 
11: Rank the solutions and find the current best. 
12: end while 

13: Output Global optimum Q g 
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there are one evaluation for the whole population in normal swarm intelligent 
methods. If we consider the number of evaluating the fitness function by these 
two evaluations, they might not be economic. 

Thus, we can make some modifications here to accelerate the CS as Al- 
gorithm 1 by decreasing the evaluation number for the fitness. 

3.2. Differential Evolution Algorithm 

Differential Evolution (DE) algorithm was proposed by Storn [79-82] . DE 
utilizes M n-dimensional vectors, X { = (xn, • • • , x in ) e S, i — 1, • • • , M, as a 
population for each iteration, called a generation, of the algorithm. For each 
vector X\ G) = (X^\x^f, • • • , X^),i — 1, 2, • • • , M, there are three main 
genetic operator acting[79-82]. 

To apply the mutation operator, firstly random choose four mutually 
different individual in the current population Xrf\Xr^\Xrf\Xrf\ri ^ 
r 2 7^ r 3 ^ f 4 ^ i) to compose a differential vector = [X^ — Xrf^] + 
[xif ~ > — Xrf\ then combines it with the current best individual xj^ s \ to get 
a perturbed vector V = (Vi, V 2: ■ • ■ , V n ) [79, 83] as below: 

V = + CF x (26) 

where CF > is a user-defined real parameter, called mutation constant, 
which controls the amplification of the difference between two individuals to 
avoid search stagnation. 

(G) 

Following the crossover phase, the crossover operator is applied on X\ . 
Then a trial vector U = (Ui, U~2, ■ ■ ■ , U n ) is generated by: 
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m 



V m , if (rand(0, 1) < CR) or (to = k), 
Xf2 , i/ (rand(0, 1) > OR) and (m^k). 



(27) 



in the current population [79], where m = 1, 2, • • • , n, the index fc e {1, 2, • • • , n} 
is randomly chosen, CR is a user-defined crossover constant [79, 83] in the 
range [0, 1]. In other words, the trial vector consists of some of the compo- 
nents of the mutant vector, and at least one of the components of a randomly 
selected individual of the population. 

Then it comes to the replacement phase. To maintain the population 

(G) 

size, we have to compare the fitness of U and X\ , then choose the better: 



3.3. Cuckoo search oriented statistically by differential evolution 

Considering the redundant evaluation for the fitness function of CS and 
the efficiency of DE, we can propose a novel cuckoo search oriented statisti- 
cally by differential evolution as following Algorithm 2. 

In each iteration of Algorithm 2, Levy flights (25) is used once for each 
location. And differential evolution operation are used with a probability 
PDe less than 0.2. In this way, the evaluations for the fitness function are 
reduced nearly 80% compared to the original Algorithm 1. 

And poe in Algorithm 2 CSODE should not be too big. Otherwise, it 
will cause the algorithm 2 be much more like a DE algorithm rather than a 
cuckoo searcher algorithm. It will be illustrated in the section simulations. 
Actually, our original idea is let CS oriented not controlled by DE. 




(28) 
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Algorithm 2 Cuckoo search oriented statistically by differential evolution 



(CSODE) 



1: 


Basic parameters' setting Objective function f(x),x = 


(xi,... 






initial population of n host nests Xi(i = 1,2, ...,n), boundaries for each 




dimension Xi = (xj,i, ...,Xi t d) and etc. 






2: 


while Termination condition is not satisfied do 






3: 


If pee < 0.2, generating candidate cuckoo population (i, i 


= 1,2, 


...,n) 




randomly from current population by equation (27). 






4: 


Updating the current cuckoo swarm and the candidate 


swarm 


with 




equation (28). 






5: 


Abandon a fraction (p a ) of worse nests. 






6: 


Build new ones at new locations via Levy flights (25). 






7: 


Keep the best solutions (or nests with quality solutions). 






8: 


Rank the solutions and find the current best. 






9: 


end while 






10: 


Output Global optimum Q g 
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4. A novel unknown parameters and orders identification method 
based on CSODE for non-commensurate fractional order chaos 
systems 

The task of this section is to find a simple but effective approach for un- 
known q and systematic parameters in equation (22) of for non-commensurate 
fractional-order chaos based on CSODE in non-Lyapunov way. 

4-1. A novel unknown parameters and orders identification method 

Now we can propose a novel approach for hyper, proper and improper 
fractional chaos systems. The pseudo-code of the proposed reconstruction is 
given below. 

Algorithm 3 A novel unknown parameters and orders identification method 

based on differential evolution algorithms for non-commensurate and hyper 

fractional order chaos systems 
1: Basic parameters' setting for Algorithm 2 . 

2: Initialize Generate the initial population. 

3: while Termination condition is not satisfied do 

4: Algorithm 2 with fitness with Eq. (22). 

5: Boundary constraints For each Xik € X i: k = 1,2, ...,D, if xn is 
beyond the boundary, it is replaced by a random number in the boundary. 
6: end while 

7: Output Global optimum XBest 



4-2. N on- commensurate and hyper fractional order chaos systems 

To test the Algorithm 3, some different well known and widely used non- 
commensurate and hyper fractional order chaos systems are choose as follow- 
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ing. To have a comparative results, these systems are taken from reference 



Example. 1. Here we discuss the non-commensurate fractional Lorenz system 
(23), so the (q u q 2 ,q 3 ) = (0.985,0.99,0.99) are taken[8]. 
Example. 2. Fractional order Arneodo's System (29) [47, 84]. 



D^x(t)=y(t); 
Dfy (t) =z(t); 

,D?z (t) = -fax (t) - fay (t) - faz (t) + fax 3 (t) . 



(29) 



when (pi, fc, Ps, fa) = (-5.5,3.5,0.8,-1.0), (gi,ga,g 3 ) = (0.97,0.97,0.96), 
initial point (—0.2,0.5,0.2), Arneodo's System (29) is chaotic. 
Example. 3. Fractional order Duffing's system (30) [47]. 



Drx{t)=y{t)- 

Dfy (t) =x{t)- x 3 (t) - ay (t) + 5 cos (wt) . 



(30) 



when (a,b,c) = (0.15,0.3,1), (q 1 ,q 2 ) = (0.9,1), initial point (0.21,0.31), 
Duffings system (30) is chaotic. 

Example. 4. Fractional order Genesio-Tesi's System (31) [47, 85]. 

' D?x(t)=y(t); 

D?y (t) =z(t); (31) 
^ Q Dfz {t) = -(3 lX (t) - (3 2 y (t) - (3 3 z (t) + fax 2 (t) . 

when (fa, fa, fa, fa) = (1.1,1.1,0.45,1.0), ( qi ,q 2 ,q 3 ) = (1,1,0.95), initial 

point (—0.1,0.5,0.2), Genesio-Tesi's System (31) is chaotic. 

Example. 5. Fractional order financial System (32) [47, 86] with the exact 
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form of the differential equation qD^z = fs(x,y,z) are unknown. 

D?x (t) = z(t)+x (t) (y (t) - a) ; 
< D?y(t) = l-by(t)-x>(t); (32) 
D?z it) = -x (t) - cz (t) . 

when (a, b, c) = (1, 0.1, 1), (q u q 2 , q 3 ) = (1, 0.95, 0.99), initial point (2, -1, 1), 
financial System (32) is chaotic. 

Example. 6. Fractional order Lu system (33) [16, 47] . 

f D?x(t)=a(y(t)-x(t)); 

< D?y(t) = -x(t)z(t) + cy(t); (33) 
^ Q Dfz(t)=x(t)y(t)-bz(t). 

when (a, 6, c) = (36, 3, 20), (q u q 2 , q 3 ) = (0.985, 0.99, 0.98), initial point (0.2, 0.5, 0.3), 
Lu system (33) is chaotic. 

Example. 7. Improper fractional order Chen system (34) [27, 47, 78]. 

D« 1 x(t) = a(y (t)-x(t)); 
< D?y(t) = (d)x(t)-x(t)z(t) + cy(t); (34) 
^ D?z(t)=x(t)y(t)-bz(t). 

And when when (a,b,c,d) = (35,3,28,-7), (gi,g 2 ,gs) = (1,1.24,1.24), ini- 
tial point (3.123,1.145,2.453), Chen system (34) is an improper chaotic 
system [10]. 

Example. 8. Fractional order Rossler System (35) [12, 47]. 

' D?x(t) = -(y(t)+z(t)); 

< D?y(t)=x(t)+ay(t); (35) 
^ Df , z{t) = b + z(t)(x(t)-c). 
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when (a,b,c) = (0.5,0.2,10), (gi,g 2 ,gs) = (0.9, 0.85, 0.95), initial point (0.5,1.5,0.1), 
Rossler System (35) is chaotic. 

Example. 9. Fractional order Chuas oscillator (36) [87]. 

D?x (t) = a(y(t)-x (t) + (x (t) - W (w) x (*)) ; 
Dfy(t)=x(t)-y(t) + z(t); 
D?z(t) = -Py(t)->yz(t); 
D?w (t) =x(t); 



(36) 



where 



W(w) 



j a : 


\w\ 




\w\ 



(37) 



when (a, /3, 7, C, a, b) = (10, 13, 0.1, 1.5, 0.3, 0.8), (q u q 2 , q 3 , q 4 ) = (0.97, 0.97, 0.97, 0.97), 
initial point (0.8,0.05,0.007,0.6), Chua's oscillator (36) is chaotic. 
Example. 10. Hyper fractional order Lorenz System (37) [88] . 

D q t 1 x(t) = a(y(t)-x(t))+w(t); 
Df 2 y (t) = cx (t) -x(t)z (t) -y(t); 
D?z(t)=x(t)y(t)-bz (t); 
^ DTw{t) = -y(t)z{t)+ 1 w{t)- 

when (a,b,c,d) = (10,8/3,28,-1), (q u q 2 , q 3 , g 4 ) = (0.96,0.96,0.96,0.96), 
initial point (0.5,0.6,1,2), Hyper fractional order Lorenz System (37) is 
chaotic. 

Example. 11. Hyper fractional order Lii System (38) [89]. 

D qi x(t) = a(y (t)-x(t))+w(t); 
D?y(t) = -x(t)z(t) + cy(t); 
D?z(t)=x(t)y(t)-bz(t); 
Q Dfw (t) =x(t)z (t) + dw (t) ; 



(38) 
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when (a, b, c, d) = (36, 3, 20, 1.3), (q u q 2 , q 3 , g 4 ) = (0.98, 0.980.98, 0.98), initial 
point (1, 1, 1, 1), Hyper fractional order Lii System (38) is chaotic. 
Example. 12. Hyper fractional order Liu System (39) [90]. 



(39) 



(40) 



Q Dfx (t) = -ax (t) + by (t) z (t) + z (t) ; 
D?y(t) = 2.5y(t)-x(t)z(t); 
Df i z(t)=x(t)y(t)-cz (t)-2w (f) ; 
Q D?w it) = -d-x{t). 

when (a,b,c,d) = (10,1,4,0.25), (q 1: q 2 , g 3 , g 4 ) = (0.9,0.9,0.9,0.9), initial 
point (2.4,2.2,0.8,0), Hyper fractional order Liu System (39) is chaotic. 
Example. 13. Hyper fractional order Chen System (40) [91]. 

Dl 1 x(t) = -a(y (t) - x (t)) + w (t) ; 
Dfy (t) = dx (t) -x(t)z (t) + cy (t) ; 
Q D?z{t)=x{t)y{t)-bz (t); 
^ D q t 4 w (t) =y(t)z (t) + rw (t) . 

when (a,b,c,d) = (35,3,12,7,0.5), (q 1: q 2 , q 3 , q*) = (0.96,0.96,0.96,0.96), 
initial point (0.5, 0.6, 1, 2), Hyper fractional order Chen System (40) is chaotic. 
Example. 14. Hyper fractional order Rossler System (41) [12]. 

Dl 1 x(t) = -(y(t) + z(t)); 
D q t 2 y(t)=x(t) + ay(t)+w(t); 

(41) 

Dfz (t) =x(t)z (t) + b; 
D q t 4 w (t) = -cz (t) + dw (t) . 

when (a,b,c,d) = (0.32,3,0.5,0.05), (91,95,93,94) = (0.95,0.950.95,0.95), 
initial point (—15.5, 9.3, —4, 18.6), Hyper fractional order Rossler System (41) 
is chaotic. 
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Example. 15. A four-wing fractional order system[92, 93] both incommen- 
surate and hyper chaotic. 

D^xi = ax\ — X2X3 + X4, 

< * 2 ' (42) 

Df x 3 = xix 2 - cx z + X1X4, 

Dfx± — —x 2 , 

when (a,b,c) = (8,40,49), (91,92,93,54) = (1,0.950.9,0.85), initial point 
(1, -2,3, 1)[92], system (42) is chaotic. 

4-3. Simulations 

For systems to be identified, the parameters of the proposed method are 
set as following. The parameters of the simulations are fixed: the size of the 
population was set equal to M — 40, generation is set to 500, the default 
values CF = 1, CR = 0.85, pde = 0.2; Table 1 give the detail setting for 
each system. 

Table 2 shows the simulation results of above fractional order chaotic sys- 
tems. And some simulations are done by single Cuckoo Search (CS) meth- 
ods. In these cases, all the other parameters for the algorithms are the same 
as for CSODE. The simulation results are listed in Table 3. 

The following figures give a illustration how the self growing evolution 
process works by DE Algorithm 3. In which, Figures 3,4 ,5,6 7,8, 9,10, 
11 show the simulation evolution results of above fractional order chaotic 
systems with optimization process of objective function's evolution and the 
parameters and orders uncertain of above fractional order chaotic systems. 

From the simulations results of above fractional order chaos system, it can 
be concluded that the proposed method is efficient. And from above figures, 
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Table 1: Detail parameters stctting for different systems 



F-0 systems Unknown Lower boundary Upper boundary Step No. of sar 



Lorenz 


(c, 7, 9i)92,93) 


1 on (i 1 ni ni ni 
0, zU, U.l, U.I, U.l, U.l 


1 r on 1 n 1 1 1 
10, OU, 1U, 1, 1, 1 


n m 
U.Ul 


1 nn 
1UU 


/\ineoQO 


//O o o a _ _ _ \ 
IPl, P2, P3, P4, 9l, 92, 93 J 


n 1 1 t: fi 1 ni (11 
— u, Z, U.l, — 1.0, U.l, U.l, U.l 


111 n 1 1 1 1 
— 0, 0, 1, — u.o, 1, 1, 1 


n nni 
u.uuo 


onn 
zuu 


Duffing 


(a, 6, c, 91,92) 


ni ni ni ni n 1 
U.l, U.l, U.l, U.l, u.o 


1, 1,2, 1, 1.5 


n nnni 
U.UUUo 


inn 
oUU 


vjrenesio- lesi 


IPl, P2, P3, P4, 9l, 92, 93 J 


1 1 ni n 1 n 1 n 1 n 1 
1, 1, U.l, U.l, u.o, U.O, U.l 


001 1111111 
z, z, 1, 1.0, 1.0, 1.0, 1 


n nni 

U.UUO 


onn 
zuu 


H 111.111 / »i o 1 


l a , °, c , 9i, 92, 93 J 


n 1 n m n 1 n 1 n 1 n 1 

U.O, U.Ul, U.O, U.O, U.l, U.l 


111 11111 1 
1.0, 1, 1.0, 1.0, 1, 1 


n nni 
u.uuo 


onn 
zuu 


T 11 


(a, 0, c, qi, q 2 , 93 J 


an n 1 k n 1 n 1 n 1 

OU, U.l, 10, U.l, U.l, U.l 


/in in 01 1 1 1 
41U, 1U, zo, 1, 1, 1 


n ni 

U.Ul 


1 nn 

1UU 


impiopei L/iieii 


(a, 0, c, a, qi, 92, 93 J 


n 1 on m ni 1 1 
OU, U.l, zu, — 1U, u.o, 1, 1 


/(l in on ni 000 
41U, 1U, ou, — U.l, Z, Z, Z 


n m 

U.Ul 


1 nn 

1UU 


Rossler 


(a b c o-i Oo cii} 


0.1,0.1,5,0.1,0.1,0.1 


1 1 15 1 1 1 

X, X, XW, X, X, X 


0.01 


100 


ChuaM 


(a, ^, 7, C, a, 9i,92,93, 94) 


5, 10, 0.1, 0.1, 0.3, 0.1, 0.1, 0.1, 0.1, 0.1 


10,20,1,2,0.3,1,1,1,1,1 


0.01 


100 


Hyper Lorenz 


(a, 6, c,d, gi,g 2 ,53,94) 


5,0.1,20,-2,0.1,0.1,0.1,0.1 


15,5,30,-0.1,1,1,1,1 


0.01 


100 


Hyper Lii 


(a, 6, c,d,q 1 ,q 2 ,q 3 ,q4) 


30,0.1,15,0.1,0.1,0.1,0.1,0.1 


40,5,25,5,1,1,1,1 


0.005 


200 


Hyper Liu 


(a, 6, c, d, gi, 92,93,94) 


5,0.5,1,0.1,0.1,0.1,0.1,0.1 


15,1.5,10,1,1,1,1,1 


0.005 


100 


Hyper Chen 


(a, 6, c, d, 91,92,93,94) 


30,0.1,10,0.1,0.1,0.1,0.1,0.1,0.1 


40,5,20,10,1,1,1,1,1 


0.005 


200 


Hyper Rossler 


(a, 6, c,d,9i,9 2 ,93,94) 


0.1,0.1,0.1,0.01,0.1,0.1,0.1,0.1 


1,5,1,1,1,1,1,1 


0.005 


200 


System (42) 


(a, 6, c, 91,92,93,94) 


5,38,45,0.5,0.5,0.5,0.5 


10,45,50,1,1,1,1 


0.005 


200 



Table 2: Simulation results for different fractional order chaos systems by CSODE 



F-0 system StD Mean Min Max Success rate a NEOF 6 



Lorenz 


2 


,2025e-04 


6 


,6862e-04 


3 


1179e- 


-06 


5 


.5594e-03 


100% 


24013 


Arneodo 


4 


,7346e-06 


1 


.5061e-05 


3 


.2217e- 


■07 


1 


.4355e-04 


100% 


24073 


Duffing 


5 


,8073e-09 


1 


,4890e-08 


1 


.3035e- 


-10 


1 


.3548e-07 


100% 


24016 


Genesio-Tesi 


1 


.3724e-04 


5. 


.1929e-04 


7. 


,6119e- 


-07 


3 


.8262e-03 


100% 


24073 


Financial 


2 


,7967e-08 


2. 


.9039e-08 


1 


,5084e- 


-09 


1 


.4728e-07 


100% 


24013 


Lu 


2. 


.7950e-04 


6 


,5471e-04 


6 


.8653e- 


-06 


5 


.9367e-03 


100% 


24038 


Improper Chen 


6 


,5650e-03 


8. 


.3188e-03 


1 


.6981e- 


-04 


4 


.3289e-02 


83% 


24074 


Rossler 


1 


.8666e-08 


2 


1368e-08 


1 


.2494e- 


-09 


1 


.2324e-07 


100% 


24013 


ChuaM 


1 


.4857e-04 


1 


.4529e-04 


1 


,9011e- 


-05 


1 


.0265e-03 


100% 


24073 


Hyper Lorenz 


5 


4741e-03 


4 


,4624e-03 


5. 


.0835e- 


-04 


2 


.4994e-02 


87% 


24041 


Hyper Lii 


1 


.2904e-02 


1 


.2573e-02 


5 


.8260e- 


-04 


5 


.5329e-02 


100% c 


24041 


Hyper Liu 


6 


,2490e-06 


4 


4206e-06 


3 


.3623e- 


-07 


2 


.3320e-05 


100% 


24101 


Hyper Chen 


5 


.4440e-02 


4 


.5710e-02 


6 


.0972e- 


-03 


2 


.7476e-01 


87% c 


24022 


Hyper Rossler 


2 


.0859e-04 


1 


5394e-04 


3 


1938e- 


-05 


7 


.9808e-04 


100% 


24041 


System (42) 


4 


.5612e-01 


1 


.8007e+00 


3 


,1952e- 


-03 


1 


.6640e+01 


56% c 


24024 



a Success means the the solution is less than f e — 2 in 100 independent simulations. 
b No. of average evaluation for objective function (NEOF) 

c Success means the the solution is less than le — 1 in 100 independent simulations. 
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Table 3: Simulation for fractional order systems by single Cuckoo Search 



system 


StD 




Mean 




Min 




Max 


Success rate a 


NEOF 6 


Lorenz 


3.8566e- 


01 


or jo 

6.8548e- 


01 


1.1083e- 


01 


2.2250 


0% 


40040 


Arneodo 


5.2993e- 


•04 


1.1322e- 


-03 


2.9726e- 


-04 


2.7731e-03 


100% 


,| A A A f\ 

40040 


Till TTI Tl CT 




D4 

U4 




D4 

U4 


7 91 87p- 




1 473flp-D3 


100% 

1UU /o 


4nn4n 


Genesio-Tesi 


4.0392e- 


03 


1.8645e- 


03 


1.0365e- 


03 


1.0422e-02 


100% 


40040 


Lii 


9.2611e- 


01 


1.8303 




6.0910e- 


01 


6.4697 


0% 


40040 


ChuaM 


8.1825e- 


02 


4.0495e- 


02 


1.5058e- 


02 


2.3472e-01 


0% 


40040 


Hyper Lorenz 


9.5917 




4.4177 




3.1812 




2.5456e+01 


0% 


40040 


Hyper Lii 


3.8130 




7.1730 




2.5631 




2.7080 


0% 


40040 


Hyper Chen 


8.8389 




3.5136 




2.6274 




2.4436e+01 


0% 


40040 


system (42) 


9.5876e- 


02 


1.4075e- 


01 


3.3331e- 


02 


5.7142e-01 


43% 


40040 



a Success means the the solution is less than le — 1 in 100 independent simulations. 



b No. of evaluation for objective function 
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Figure 4: Evolution process for fractional order Arneodo system 




Figure 5: Evolution process for fractional order Financial system 
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Figure 6: Evolution process for fractional order Lii system 
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Figure 8: Evolution process for fractional order hyper Lii system 




Figure 9: Evolution process for fractional order hyper Liu system 
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Figure 11: Evolution process for fractional order four wing system (42) 
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it can be concluded that the estimated systems are self growing under the 
genetic operations of the proposed methods. 

To test the performance of the proposed method Algorithms 3 , some 
more simulations are done to the four-wing incommensurate hyper fractional 
order chaotic system (42) in following cases A,B,C,D. In these cases, each 
with only one condition is changed according to the original setting for system 
(42). The other parameters for the algorithms are the same as for CSODE. 
The simulation results are listed in Table 4. 

• Case A. Enhancing the defined intervals of the unknown parameters 
and orders to [0, 10] x [30, 45] x [40, 50] x [0.1, 1] x [0.1, 1] x [0.1, 1] x [0.1, 1] 

• Case B. Reducing the number of samples for computing system (42) 
from 200 to 100. 

• Case C. Increasing the iteration numbers of Algorithms 3 from 500 to 
800. 

• Case D. Changing the population size of Algorithms 3 from 40 to 80. 

Figure 12 show the coresspondent simulation results for system (42). 

From results of the Table 2,3, 4 and Figure 12, we can conclude that 
minimizing the number of samples for computing the system (42) as case 
B, enhancing the iteration numbers as case C, the population size of Al- 
gorithms 3 as case D, will make the Algorithms 3 much more efficient and 
achieve a much more higher precision. However if the defined intervals of the 
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Table 4: Simulation results for system (42) 



system (42) 


StD 


Mean 


Min 


Max 


Success rate 


NEOF 6 


Case A. 


2.3271e+02 


2.2944e+02 


3.4062e-02 


5.2168e+02 


31% a 


24000 


Case A. c 


3.5361e+02 


4.3992e+02 


2.9705e-04 


7.2477e+02 


39% a 


24040 


Case B. 


3.3699e-02 


1.6530e-01 


1.5989e-05 


8.3946e-01 


96% d 


23800 


Case C. 


6.2779e-03 


2.4339e-02 


1.4415e-05 


1.7882e-01 


98% d 


39360 


Case D. 


7.5365e-01 


2.3014e+00 


1.1463e-02 


1.6874e+01 


41% d 


48320 



a Success means the the solution is less than 1 in 100 independent simulations. 
b No. of evaluation for objective function 
c By single DE with 600 generations 

d Success means the the solution is less than le — 1 in 100 independent simulations. 







Figure 12: Simulation results for system (42) 
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unknown parameters of system (42) are enhanced, then the results will go to 
the opposite way. That is the success rate is from 90% to 20% as case A. 

If the No. of evaluation for objective function is considered, it is that 
minimizing the number of samples for computing the system (42) as case B 
is the best way to achieve higher efficiency and precision. 

And according to Table 2,3, the Algorithms 3 based on CSODE is much 
more better than single cuckoo search. 

5. Conclusions 

The novel Gao's mathematical model in Section 2 is not only for fractional 
chaos parameters identification but also for reconstruction in non-Lyapunov 
way with three sub-models A, B and C. 

The put method based on CSODE consists of numerical optimization 
problem with unknown fractional order differential equations to identify the 
chaotic systems. Simulation results demonstrate the effectiveness and ef- 
ficiency of the proposed methods with the Gao's mathematical model in 
Section 2. This is a novel Non-Lyaponov way for fractional order chaos' un- 
known parameters and orders. The proposed the method solve a the question 
that the unknown fractional order q are not resolved in reference [36]. 

Here we have to say that this work is only about the estimation of un- 
known parameters and orders with the objective function (22) for for non- 
commensurate and hyper fractional order chaos systems in non-Lyapunov 
way. It can be concluded that CSODE in Algorithms 3 can be change to 
other artificial intelligence methods easily. And from the Table 2,3, 4 , we 
can conclude that Algorithms 3 with CSODE is better than CS. And the the 
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system (42), CSODE is better than DE for the bigger scale of the unknown 
parameters and orders. 

In the future, there are three interesting problems to be studied. 

• Neither the fractional orders nor some fractional order equations are 
unknown. That is, the objective function is chosen as (12) in the novel 
mathematic model in Section 2. A simple way for this might be the 
approaches combining the fractional orders and fractional order equa- 
tions together, that might be both the estimation methods as artificial 
intelligent methods and the reconstruction methods as in Reference [36] 
together in some degree. 

• Time-delays and systematic parameters are unknown for fractional 
time-delay chaos systems. The objective function will be selected as 
the objective function (22) for the time-delay fractional chaotic systems 
as in Gao's sub-model C , which have special characteristics. 

• Cases with noises. Normally, the white noise will be added to the Gao's 
three sub-models. The similar ideas discussed in Section 2 will be used. 

In conclusion, it has to be stated that proposed Algorithms 3 for fractional 
order chaos systems' identification in a non-Lyapunov way is a promising 
direction. 
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